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O ; ABSTRACT 

, The galaxy population at z « 6 has been the subject of intense study in recent 

■ years, culminating in the Hubble Ultra Deep Field (HUDF) - the deepest imaging 

\ survey yet. A large number of high redshift galaxy candidates have been identified 

within the HUDF, but until now analysis of their properties has been hampered by 
I ' the difficulty of obtaining spectroscopic redshifts for these faint galaxies. Our "Gemini 

Lyman- Alpha at Reionisation Era" (GLARE) project has been designed to undertake 
. spectroscopic follow up of faint (z! < 28.5) z'-drop galaxies at z ss 6 in the HUDF. 

In a previous paper we presented preliminary results from the first 7.5 hours of data 
. . , from GLARE. In this paper we detail the complete survey. We have now obtained 36 

hours of spectroscopy on a single GMOS slitmask from Gemini-South, with a spectral 
| resolution of A/AAfwhm ~ 1000. We identify five strong Lyman-a emitters at z > 5.5, 

. and a further nine possible line emitters with detections at lower significance. We also 

' place tight constraints on the equivalent width of Lyman-a emission for a further ten i - 

drop galaxies and examine the equivalent width distribution of this faint spectroscopic 
sample of z « 6 galaxies. We find that the fraction of galaxies with little or no emission 
is similar to that at z ~ 3, but that the z « 6 population has a tail of sources with high 
rest frame equivalent widths. Possible explanations for this effect include a tendency 
towards stronger line emission in faint sources, which may arise from extreme youth or 
low metallicity in the Lyman-break population at high redshift, or possibly a top-heavy 
initial mass function. 
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1 INTRODUCTION 

The Hubble Ultra Deep Field (HUDF. lBeckwith et alj|2006h 
opened a new window on the early universe. The exception- 
ally deep, multiwavelength data provided opportunities to 
study colour-selected samples of high redshift galaxy candi- 
dates with modest luminosities more typical of the general 
galaxy population. Previously, only the most luminous 'tip 
of the iceberg' had been accessible. By pushing to redshifts 
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around 6 with the z'-ban d filter and the i'-dr op Lyman break 
selection technique (e.g. iBunker et al]|2004l ). the HUDF ex- 
plored the end of the reionization epoch sig naled by the 
Gunn -Peterson HI absorption trough in QSOs dBecker et al.l 
l200lfi . 

Nonetheless, the use of the Lyman break colour se- 
lection criterion to isolate star-forming galaxies at a spe- 
cific redshift, first developed to study z w 3 galaxies 
jSteidel et all Il995h and extended to z ~ 6 candidates 
in analysis of the Great Observatorie s Origins Deep Sur- 
vey (GOODS) and HUDF fields (e.g. IStanwav et al.ll2003l ; 
iBouwens et al.l [2006t ) . presents challenges. Before the prop- 
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erties of the population can be meaningfully discussed, the 
selection function itself must be understood. Estimates must 
be obtained of the contaminant fraction and the redshift dis- 
tribution of i'-drop galaxies. 

Many colour-selected galaxies are significantly fainter 
than the conventional spectroscopic limit of today's large 
telescopes. The HUDF reaches limits of z' AB = 28.5 (10 a) - 
a depth that, until JWST becomes available, requires unrea- 
sonable exposure times to obtain a high signal-to-noise ratio 
(S/N) spectrum for the continuum. However, line emission 
(e.g. Lyman-a) may be within reach of ultra-deep moderate- 
dispersion spectroscopy on 8-10 m telescopes even for the 
faintest galaxies, provided the equivalent width of the line is 
large enough. If the properties of the star-forming popula- 
tion at z w 6 are similar to th ose of the Lyman Br eak Galaxy 
(LBG) population at z « 3 (|Steidel et alj|2003j ). then half 
the sources are expected to show the resonant Lyman-a 
transition in emission. Detection of a Lyman-a emission line 
fixes the redshift of a source, while detection or constraints 
on nearby high ionisation emission lines can quantify contri- 
bution by an AGN. Furthermore, spectroscopy can enable 
the identification of lower-redshift contaminants in the sam- 
ple which will not emit Lyman-a photons at these wave- 
lengths, but may show other emission lines. When a pop- 
ulation of Lyman-a emitters is studied, the distribution of 
equivalent widths is sensitive to the stellar initial mass func- 
tion (since the Lyman-a transition is excited by emission 
from hot, massive, short-lived stars), and also to the frac- 
tion of neutral gas in the inte rgalactic medium (IGM) an d 
the emergent photon fraction (|Malhotra fc R hoads 2004b). 

In this paper we present results from the Gemini 
Lyman-Alpha at Reionisation Era (GLARE) project. This 
program used the 8-m Gemini South telescope to obtain 
extremely deep spectroscopy on a single slitmask, centered 
on the HUDF (Figure [T]). By obtaining extremely long expo- 
sures using a telescope with a large collecting area, we aimed 
to study continuum-selected galaxies fainter than those tar- 
geted by any other survey, and to quantify the line emission 
properties of our i'-drop sample in th e HUDF. We presented 
initial results from this programme in lStanwav et al] (|2004bl . 
hereafter Paper I) , which was based on 7.5 hours of on-source 
exposure taken at low spectral resolution (A/ A A m 500). In 
this paper we present an analysis of the twenty-four i'-drop 
selected 5.6 < z < 7.0 candidates targeted for 36 hours of 
spectroscopy at higher resolving power (A/AA « 1200). In 
section [2] we describe the GLARE program, and in section 
[3] we present the results of our observing campaign. In sec- 
tion 2] we analyse the equivalent width distribution of the 
GLARE line emitters, and in section[5]we discuss the impli- 
cations of this distribution in the context of other work in 
this field. Finally in section [5] we present our conclusions. 

We adopt the following cosmology: a flat Universe with 
n A = 0.7, Q M = 0.3 and H = 70/i7okms" 1 Mpc" 1 . All 
magnitudes are quoted in the AB system (Oke & Gunn 
1983). 



2 THE GLARE PROJECT 
2.1 Candidate Selection 

The objectives of the GLARE project were to determine the 
redshift and line emission properties of our targets, reach- 
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Figure 1. The geometry of the UDF and GOODS-S fields (taken 
at two position angles), and the layout of the GLARE mask. 
GLARE encompasses the entire UDF field, and some parts of the 
GOODS-S field 



ing low equivalent width limits on Lyman-a emission for 
a u niformly-selected sample of i'-drop Lyman-break galax- 
ies llStanwav et al]l20o3 : iBouwens et al.ll2004l : IBunker et aD 
I2004T ). The i'-drop colour selection is sensitive to galaxies 
at 5.6 < z < 7.0, although the selection sensitivity falls 
off rapidly above z ps 6.5 (see Stanway et al. 2003). With 
this data we aim to probe galaxies at the end of the epoch 
of reionisation, and compare their properties to comparable 
galaxy samples at lower redshift, in order to explore possible 
evolution in the nature of the star-forming galaxy popula- 
tion. 

Initial results obtaine d from this survey were reported 
in IStanwav et af] (|2004bl ). which presented 7.5 hour spec- 
tra of three Lyman-a line emitters observed on the 2003 
GLARE slitmask. Selection of targets for our 2003 observa- 
tions was made possible by the early release of a list of red 
sources in the HUDF, based on observations to one third 
of its final depth, and supplemented by i'-drop galaxies se- 
lected from the shallower, but wider field, GOODS survey 
of the same region (Stanway 2004). 

Between observations of our 2003 GLARE mask, and 
the start of semester 2004B, the full HUDF data were made 
public and the GLARE slitmask was redesigned accord- 
ingly. Sources with (i' — z')ab 1.3 and z' AB ^ 28.5 in 
the HUDF imaging (i.e. the catalog of IBunker et al.ll2004T l 
were assigned the highest priority for spectroscopic follow- 
up and fifteen such sources were placed on the slitmask. 
The agreement is excellent between t he HUDF i'-drop dis- 
covery catalog of IBunker et al] (|2004h . which we use as our 
GLAR E source list, and the recent data paper by the HUDF 
team l|Beckwith et al] 120061 ) . All the robust i'-drops tar- 
getted on the G LARE slitmask were are also identified by 
iBeckwith et al] 

In order to maximise use of the mask, additional ob- 
jects were targetted. Seven slits were placed on candidates 
(including the two emitters identified in Paper I, GLARE 
3001 & 3011) from a brighter selection with (i! — z')ab ^ 1-3 
and z' AB ^ 27.5, selected from the GOODS imaging (Stan- 
way 2004), primarily at the edges of the slitmask and ly- 
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ID 



No of Slits Description 



1XXX 15 HUDF, i' - z' ^ 1.3, «' ^ 28.5 

2 HUDF, other red sources 

2XXX 5 Alignment Stars 

3XXX 7 GOODS i' - z' ^ 1.3, z' < 27.5 

4XXX 2 GOODS v - i' selected, 2 = 5 cands 

5XXX 2 Known, low z, [Oil] emitters 

6XXX 6 Blanks sky slits 

7XXX 3 GOODS marginal i'-drop candidates 



12 



Total Number of Slits on Mask 



Table 1. The composition of the 2004 GLARE slitmask. 1XXX 
indicates a GLARE identifier number in the range 1000-1999, etc. 



ing outside the HUDF region. A further three slits were 
placed on candidates with colours lying marginally out- 
side of our GOODS selection criteria, either in colour (i.e. 
1.0 < i' — z' < 1.3) or magnitude (i.e. z > 27.5), or in noisy 
regions of the GOODS images, and two slits on sources from 
the original UDF early release list of red sources that did 
not meet our final criteria. Two slits were placed on z ps 5 
u-band dropout candidates, and two on known low redshift 
[Oil] 3727 A emitters , previously observed by VLT/FORS2 
l|Vanzella et al.l [20061 ) and used as a check on flux calibra- 
tion. Finally, five slits were placed on alignment stars to en- 
sure accurate positioning of the mask, and six slits were used 
to conduct a blank sky survey for serendipitous sources. The 
composition of the 2004 GLARE slitmask is summarised in 
table[T] 

Both the z' AB = 27.5 cut applied in the GOODS data 
and the z' AB = 28.5 cut applied in the HUDF correspond 
to a signal to noise of approximately 8. We chose to work 
well above the detection limit in order to have confidence 
in the reality and nature of our candidate sources (many of 
which are detected only in this one band). In the event of 
non-detection in the i'-band we use the limiting magnitudes 
i' AB = 29.15 (GOODS) and i' AB = 30.4 (HUDF) correspond- 
ing to 3 a variation in the sky background, as measured on 
the images. All sources were required to be undetected in 
the available b (F435W) imaging but faint detections in the 
deep v (F606W) filter (which lies above the Lyman limit 
at z « 6, and which is present in several spectroscopically 
confirmed z > 5.6 galaxies) were permitted. Completely un- 
resolved sources were omitted from the selection although at 
faint magnitudes, the dividing line between unresolved and 
slightly-resolved sources becomes blurred. 

Of the line emitter candidates presented in section [31 
only GLARE 1042 and GLARE 1040 lie within the NIC- 
MOS HUDF field (|Thompson et al.ll2005h . Both have near- 
infrared colours consistent with a high redshift interpreta- 
tion, as i ndeed did all the i'-d rop sources in this field dis- 
cussed bv lStanwav et all l|2005h . 



2.2 Observations and Data Analysis 

The 2004 GLARE slitmask was observed in semester 2004B 
using the GMOS spectrograph on Gemini-South (Hook et 
al. 2003) . This mask was observed at higher spectral resolu- 
tion than the 2003 GLARE mask, using the GMOS-S R600 
grating rather than the R150, blazed at an angle of 48.5° 
giving a central wavelength of ~ 8700A. As the CCD over- 



sampled the typical seeing and spectral resolution, the image 
was binned at 2 x 2 pixels so as to reduce the readout noise 
and improve the S/N. After this binning, the spectral scale 
was 0.94 A/pixel, and 0.146 arcsec/pixel spatially on the de- 
tector. The slit width was 1.0 arcsec, which produced a spec- 
tral resolution of 6.5 A FWHM for objects which fill the slit. 
Both mask and slits were oriented due north. The higher 
spectral resolution of the 2004 GLARE mask (A/ A A « 1200 
compared with 500 for 2003 GLARE) decreases the fraction 
of the wavelength range adversely affected by OH sky lines, 
and also enables us to better study the profiles of emis- 
sion lines from the targets. The OG515 filter was used to 
suppress second order light from shorter wavelengths. In or- 
der to fill CCD chip gaps and ensure full wavelength cov- 
erage in the range A » 7000 - 10000A three different cen- 
tral wavelength settings were observed (8580A, 8700A and 
8820 A). The shortest wavelength reached was 6420 A, while 
the longest wavelenth surveyed was 10950A. 

Wavelengths were calibrated from the night sky lines in 
each slit, leading to a solution with an rms error of approx- 
imately 0.3A. Fluxes were calibrated from the broadband 
magnitudes of the alignment stars on the mask, and checked 
against both line emission of known [Oil] emitters also ob- 
served in VLT/FORS2 spectroscopy that had been placed 
on our mask for verification purposes, and existing spec- 
trosc opy for the known z = 5. 83 Lyman-q emitter GLARE 
1042 jStanwav et al.ll2004bl lal; lDickinson et alH2004T ). We es- 
timate a 20% error on the flux calibration, associated with 
slit losses and centroiding uncertainty in the narrow slits. 

To optimise the subtraction of night sky emis- 
sion lines (which occupy a large fraction of the spec- 
trum at > 8000A) we used the instrument in 'nod 
fc shuffle' mode fclazebrook fc Bland-Hawthorn! l200ll ; 
I Abraham et all 120041 '). Each exposure was 30 minutes long, 
nodding every 60 seconds. Hence we are able to suppress sky 
emission that varies on timescales longer than one minute. 
The total exposure time on this mask was 36 hours. 

The reduction of nod & shuffle data involves the sub- 
traction of positive and negative spectra, observed in al- 
ternate 60 second exposures and offset spatially by the 
telescope nodding. We employed slits 2.47 arcseconds in 
length, nodding by 1.25 arcseconds between sub-exposures. 
Our queue scheduled observations were constrained by the 
requirement that the seeing was less than 0.5 arcseconds 
FWHM, and the nod distance set such that the signals were 
separated by at least twice the seeing disk. Hence the charac- 
teristic signal of an emission line comprises a 'dipole' signal 
of positive and negative emission at the same wavelength, 
spatially offset by 1.25 arcseconds. In visually identifying 
candidate line emitters, we looked for this dipole signature 
with positive and negative channels of comparable strengths; 
this requirement effectively increases the sensitivity of the 
spectroscopy beyond the formal limit, since random back- 
ground fluctuations are unlikely to produce simultaneous 
signals in both positive and negative channels. 

Dipole signals lying under strong emission sky lines are 
treated with caution; as well as having more Poisson noise, 
sky line variability and charge diffusion at the red end of the 
spectrum may lead to a spurious signal. 

A number of charge traps and CCD artifacts were 
also masked when the individual exposures were combined. 
These charge traps affect some regions of the CCD more 
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severely than others, and so the noise properties vary from 
slitlet to slitlet, and also with wavelength. However, the 
1 sigma rms pixel-to-pixel variation in the background at 
8500A was 1.4 x 10" 19 erg cm -2 s _1 A -1 for the 2 x 2 binned 
pixels, measured between skylines on the nod & shuffle 
background-subtracted spectra. Hence, our sensitivity to 
an emission line extending over 500 km s" 1 x 1 arcsec is 
2.5 x 10 -18 erg cm -2 s _1 (3 a combining both nod positions, 
or 2a per nod channel). 



3 LINE EMITTERS IN THE 2004 GLARE 
MASK 

In the 36 hour exposure of the 2004 GLARE mask, we iden- 
tify five strong Lyman-a emission line sou rces (including 
the th ree candidates tentatively proposed in IStanwav et all 
l2004bl ). We identify a further four sources which have lower 
significance detections, but which are considered possibly to 
be real emission lines due to their dipole natures and sep- 
aration from possible sky line residuals. Finally we identify 
five sources with tentative emission line detections that are 
considered unlikely to be real. 



3.1 Strong Emission Lines 

Table [5] lists the properties of the five z ~ 6 sources 
with strong Lyman-a line emission. Figure [2] presents the 
two dimensional spectra obtained for these sources, and 
the summed flux from positive and negative spectral chan- 
nels. We also placed two known lower-redshift galaxies on 
the 2004 GLARE mask, placing [Oil] within our spec- 
tral range, as a check on our sensitivity and calibration. 
These are galaxies GDS J033235.79-274734.7 (z = 1.223) 
and GDS J033229.63-27 4511.3 (z = 1 33) fr om the ESO 
VLT/FORS2 survey of IVanzella et all i|2005D . We detect 
the [Oil] emission at A = 8285 A & 7577 A, with line fluxes 
1.53 & 1.52 x 10 -17 ergcm~ 2 s" 1 , respectively. 

For the Lyman-a detections, GLARE sources 1042 and 
3001 were previously identified as line emitters i n Pape r 
I from the 2003 GLARE mask, IStanwav et"aH (|2004bl ). 
GLARE 3011 was tentatively identified as a line emitter, 
and this identification is now strongly confirmed. GLARE 
1054 and 1008 are presented for the fi rst time here. 

Since its discovery (SBM03#1 in IStanwav et~al]|2003r i 
GLARE 1 042 at z = 5.83 h as been spectroscopica lly 
confirmed (|Stanwav et alj l2004al lbt iDickinson et all 120041 - 
SiD002) and extensively studied, including in t he infrared 
with Spitzer (|Evles et al.ll2005l ; Fan et alj|200l - #lab). 
GLARE 1054 was identified from the i'-drop selection in the 
HUDF, GLARE 3001 and 3011 from the somewhat brighter 
GOODS selection. GLARE 1008 is a source selected from 
the initial early data release list of red sources in the HUDF 
survey. It lies outside the final overlap region of the HUDF 
which has been used for catalogue construction and analysis, 
but within the noisy outer regions of the HUDF mosaic. It is 
technically below the detection limit of the GOODS survey, 
although it is faintly detected in the GOODS z'-band. It is 
faintly detected in the z'-band in the shallower edges of the 
HUDF. 

Given the higher resolution and more sensitive flux limit 



of this spectroscopy, we can rule out a low redshift interpre- 
tation for these sources and confirm them as z m 6 Lyman-a 
emitters. The spectral resolution of the GMOS configuration 
used in our 2004 GLARE mask is sufficient to resolve the 
[Oil] A rost = 3727, 3729A doublet (and does so in the case of 
the two known [Oil] emitters on our 2004 GLARE mask), 
and we would expect to identify weaker emission lines such 
as [Nil], [SII] within our observed redshift range if the de- 
tected line was Ha (A re st=6563A). The other strong opti- 
cal emission lines (H/34861, [OIII] 4959,5007) would always 
yield other strong lines within our spectral range. Sources 
at these low redshifts are also unlikely to satisfy the colour 
cut criterion used for target selection. 

All five detected strong Lyman-a emission lines also 
show significant asymmetry as expected for high redshift 
emitters (in which the blue wing of the line is significantly 
self absorbed, and the red wing broadened by re-emission). 
This phenomenon is well known at z ~ 3 and is believed 
to a rise in outflows from t he actively star-forming galax- 
ies (|Adelberger et all 120031 ). Similar asymmetry has now 
been observed in all z > 5.5 galaxies confirmed to date by 
spectroscopy that resolves the Lyman-a emission line (e.g. 
iBunker et alll2003l ). suggesting that similar outflows are pro- 
duced at higher redshift galaxies. 

3.2 Possible Emission Lines 

We identify a further four sources (listed in table |3J) as pos- 
sible line emitters. In most cases, the detection in each chan- 
nel is of low significance, but the coincidence of positive and 
negative signals suggests that the emission lines are real. 
Alternately some candidate lines may lie on top of skylines, 
leading to the suspicion that these represent skyline resid- 
uals. These sources are illustrated in figure [3] All but one 
of the candidate emission lines lie shortward of the lower 
redshift limit selected by the i'-drop technique. Sources at 
these redshifts would be expected in a w-drop rather than 
i'-drop selection. It is possible for sources with large errors 
on their i' — z' colour to scatter into the i'-drop selection, 
although it is unlikely that photometric scatter alone could 
explain this discrepancy. 

Three of these sources were identified in the HUDF i! - 
drop sample, one (GLARE 3000) from the GOODS sam- 
ple. GLARE 3 00 wa s first identified as an i'-drop source in 
IStanwav et~ai1 (|2003l . , SBM03#05) but considered likely to 
be an M or L class Galactic dwarf star, as it is unresolved 
in HST imaging. The candidate emission line in this source 
lies beside a stro ng sky line res i dual. FORS2 spectroscopy 
of this source bv lVanzella et a l. (2004) also interpreted the 
spectrum obtained as that of a Galactic star. The stellar 
identification is almost certainly correct, given the broad- 
band colours and unresolved half light radius of this source, 
suggesting that possible emission lines at this significance 
should be considered highly suspect. 

There are two possible emission lines of similar strength 
in GLARE 1067 (at 7037A & 7099 A), rendering it unlikely 
that this source lies at high redshift. The observed line sep- 
aration may be consistent with [OIII] emission (A rcst = 
4959, 5007A) at z = 0.418, although a galaxy at this red- 
shift is predicted to have colours that are significantly bluer 
(i.e. i' — z' < 0.5 and detected in the b band, as opposed 
to the observed i' — z' = 1.4 ± 0.2 and no b detection). In 
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Figure 2. One and two dimensional spectra of the strong emission line candidates. The upper panel (a) shows the two dimensional 
spectrum, the central panel (b) a spectrum smoothed over three pixels, and the lower panel (c) the one dimensional combined spectrum 
extracted from both positive and negative channels. The source ID is shown on the right of panel (c) in each case. The 1 sigma standard 
deviation in the background, smoothed over 15A and including the Poisson noise on the emission line flux is indicated with a dotted 
line. The arrow indicates the wavelength of the putative emission line. 



ID 


RA & Declination 
J(2000) 


Alternate 
ID 


z 


z 1 

AB 


i< -a! 


Line Flux 
erg cm -2 s _1 


w rcst 

Lya 
A 


^Lya 

10 42 ergs" 1 


1042t 


03 32 40.01 -27 48 14.9 


20104 


5.83 


25.35±0.02 


1.64±0.04 


1.58 xlO" 17 


22.9 


5.90 


1054 


03 32 33.20 -27 46 43.3 


42414 


5.93 


27.65±0.07 


1.45±0.17 


0.68 xl0~ 17 


120 


2.63 


1008 


03 32 47.97 -27 47 05.1 




6.13 


28.51±0.18 


>1.35 


0.43 xl0~ 17 


159 


1.80 


3001 


03 32 46.04 -27 49 29.7 




5.79 


26.69±0.06 


1.67±0.20 


0.77 xlO" 17 


44.1 


2.83 


3011 


03 32 43.18 -27 45 17.6 




5.93 


27.47±0.12 


1.86±0.50 


1.13 xlO" 17 


242 


4.39 



t GLARE 1042 has a panoply of alternate names. It is SBM03#3 in lStanwav et alJ d2003h. 2 0104 in l|Bunker et al.ll2004l ) and is 

SiD0002, the z = 5.83 line emitter of lDickinson et all fcOoH l. 

Table 2. Strong line emitters on the 2004 GLARE mask. Errors on line fluxes and equivalent widths are approximately 20%. Equivalent 
widths are calculated from broadband magnitudes, accounting for line contamination a nd Lyman-o forest blanketing. 2 a limits on 
magnitudes are given where appropriate. Alternate ID is taken from Bunker et al ] l|2004h . Note that GLARE 1042 lies in the chip gap 
of one of the three wavelength settings observed, and hence in a region with slightly lower signal to noise. 



addition, the redward line of this [OIII] doublet is expected 
to be three times stronger than the blueward line, while the 
observed emission peaks are of comparable strength. An al- 
ternate explanation might be that two galaxies, separated by 
3000 kms -1 lie within the slit. While this source has a close 
neighbour (at a.12000 =03:32:35.8, <5.j20oo =-27:48:49, with 
a separation of < 1 arcsecond), the neighbouring galaxy is 
blue in colour (i 1 — z' = 0.2), and likely lies at significantly 
lower redshift. Therefore the most probable interpretation 
of these lines (if they are real) is that they represent [OIII] 
emission arising not in the targeted galaxy but rather in the 
neighbouring source. 

GLARE 1040 and GLARE 1086 are the most convinc- 
ing candidates in this category. GLARE 1040, is an isolated 
source that clearly drops between the i' and z'-bands, and 
is undetected in the HUDF -u-band. The candidate emission 
line lies at z = 5.2, just shortward of the i'-drop selection, al- 
though it is possible for a faint galaxy such as this to scatter 
upwards into the selection window. 

The candidate line in the final source, GLARE 1086, lies 



firmly in the redshift range selected by the i'-drop technique, 
and may well be a high redshift emission line, although the 
detection is too weak to rule out a low redshift [Oil] expla- 
nation, or to detect line asymmetry. 



Finally, we identify signals in a further five sources as 
'unlikely' emission lines. These are shown in table U and 
figure[4] Although the dipole signal from the first two sources 
appears strong, they lie on bright skylines and so may be 
partly or wholly due to sky subtraction residuals. They may 
also arrise as a result of charge diffusion from the adjacent 
slit. Both candidates lie at A > 10000A (not far from the 
10500A Si bandgap) where the diffusion is strong. The signal 
in the remaining candidates is weak. Slit GLARE 6050 was a 
blank sky slit, and there is nothing visible in any waveband 
at the depth of the GOODS imaging of this region. The 
possible serendipitous line emitter is offset from the centre 
of the slit by approximately (f.'2 to the north. 
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Figure 3. One dimensional combined spectrum extracted from both channels of the image, and two dimensional spectra, of the possible 
emission line candidates. As in figure [2] 



Slit 


RA &; Declination 
J(2000) 


Alternate 
ID 


z 


AB 


i'-z' 


Line Flux 
erg cm -2 s — 1 


wrest 
Lya 

A 


10 42 ergs" 1 


1040 


03 32 38.28 -27 47 51.3 


24458 


5.21? 


27.51±0.07 


1.60±0.17 


0.17X10" 17 


20.2 


0.49 


1067 


03 32 35.83 -27 48 48.9 


14210 


4.84?t 


28.08±0.10 


1.43±0.24 


0.30xl0~ 17 


81.5 


0.72 








4.78?+ 






0.31X10" 17 


88.0 


0.73 


1086 


03 32 30.14 -27 47 28.4 


30359 


6.10? 


28.13±0.11 


1.46±0.25 


0.37xl0~ 17 


68.1 


1.53 


3000* 


03 32 38.80 -27 49 53.7 




5.08? 


25.65±0.03 


>3.50 


O.llxlO" 17 


2.24 


0.30 



t GLARE 1067 has two emission lines, which are consistent with [OIII] at z = 0.418; if these are instead Lyman-a, both are at z < 5. 

* GLARE 3000 is most likely a low-mass Galactic star. 

Table 3. Possible line emitters on the 2004 GLARE mask. As in table[2] The redshift assumes that the emission detected is the Lyman-ct 
line at A rcst = 1215. 7A. 
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Figure 4. One and two dimensional spectra of the 'unlikely' emission line candidates presented in table [4] As in figure [2] 
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Slit RA & Declination Alternate 
J (2000) ID 



1034 

1030 

3015 

6050* 

8005 



03 32 39.454 -27 45 43.42 
03 32 41.184 -27 49 14.81 
03 32 27.910 -27 49 41.98 
03 32 37.450 -27 49 47.60 
03 32 34.392 -27 45 33.01 



52086 7.24? 27.97±0.09 

10188 7.46? 27.10±0.05 

5.67? 27.30±0.10 

7.05? 

4.69? 27.11±0.16 



>2.43 
2.04±0.16 
>1.85 

1.15±0.36 



6050 is a 'blank sky' slitlet 



Txrrcst 
Lya 
A 



10.6 

4.4 

16.7 

>104 

66.6 



Table 4. The five emission line candidates classified as 'unlikely' to be real. A limit has been placed on the possible serendipitous emitter 
GLARE 6050 based on non-detection in the z' band to our imaging depth of z' Ag = 28.5. Column headings as in table [2] 



3.3 Sources with No Line Emission 

The remaining 21 science targets on the slitmask showed no 
evidence for line emission in the wavelength range observed, 
Aobs ~ 7800 - 10000A (corresponding to 5.4 < z < 7.2, with 
some slit to slit variation depending on slit location). 

Of these, 11 formed part of our high priority i'-drop 
selection, with marginal targets (5 sources) and sky slits (5 
slits) constituting the remainder. We discuss the implica- 
tions of these non-detections of line emission below. 

It is possible that, despite optimising our experimen- 
tal setup for sky line subtraction, we are missing flux from 
emission lines that fall in regions dense with skylines. The 
noise in such areas is greater than the slit average, making 
identification of line candidates more difficult. Some 35% of 
the wavelength coverage of this spectroscopy lies under sky 
lines (defined as regions where sky emission exceeds 150% of 
the inter-line sky continuum level) which leave Poisson noise 
and sky brightness fluctuation residuals of varying strength, 
although only 10% lies under strong skylines (sky line emis- 
sion > 5 times sky continuum). The noise under the resid- 
uals of bright sky lines is a factor of 2.5 greater than that 
between lines, and might also be prone to systematics in 
the sky subtraction leading to potentially spurious emission 
lines. 

Only half of Lym an Break Galaxies at z ~ 3 show 
Lyman- a in emission (jShaplev et al.l [20031) . and one quar- 
ter have W^yl >20A. Given that we reach a limiting rest- 
frame equivalent width of ^6A for the majority of our tar- 
gets, we might naively expect to observe emission in half our 
slits, less ten percent lost to strong sky line residuals. How- 
ever, the fraction of Lyman Break Galaxies with emission 
at 2 w 6 is still unknown, and it is likely that a large frac- 
tion of z ps 6 i'-drop galaxies will not be spectroscopically 
confirmed until the continuum level is reached with either 
longer exposures or more sensitive telescopes, and interstel- 
lar absorption lin es can be used for redshift determination 
(see, for example, iDow-Hygelund et alJliooBl for rest-frame 
UV continuum spectroscopy of a bright z = 5.5 galaxy). 

Finally in the case of the marginal candidates used to 
fill the slitmask, it is possible that the true redshift of the 
sources lies above z = 4.5 (set by the dual requirements of 
non-detection in the b band, and the i' — z' colours of these 
sources), but below the limit of our spectroscopy. Two of 
the slits without detections were placed on v-dmp galaxies 
expected to lie in this range, while a further four slits were 
placed on sources slightly too blue to meet our strict se- 
lection criteria. It seems unlikely that photometric scatter 
could place these sources in our detectable redshift range. 



We will discuss the equivalent width limits on our i! - 
drop galaxies in Section 2] 



3.4 Agreement with Other Spectroscopic Surveys 

A subset of sources in our sample have also been observed 
spectroscopically as part of the concerted campaign of follow 
up observations to the GOODS survey. 

Sources in the GOODS-S field have been targeted for 
8m spectroscopy by FORS 2 on the VLT (Vanze lla et al, 
2006), by DEIMOS on Keck l|Stanwav et al.ll2004al . , Bunker 
et al 2006) and by GMOS on Gemini (Paper I and this work). 
In addition, this field was surveyed with the HST/AC S 
grism as part of the GRAPES survey dMalhotra et alll2005l ). 
This slitless spectroscopy has been obtained to varying 
depths, and at different spectral resolutions. In particular, 
the GRAPES grism spectroscopy is of too low a resolu- 
tion to detect all but the highest equivalent width Lyman-Q 
emission lines adjacent to the Lyman-a forest continuum 
break (compare, for example, the DEIMOS spectrum of 
SBM03#1/GLARE1042 in Stanway et al. 2004 showing the 
strong Lyman-a line, with the HST/ACS spectrum of the 
same object, SiD002, in fig. 2a of Dickinson et al. 2004, which 
shows a continuum break but the line is washed out). How- 
ever, GRAPES does provide a redshift estimate by localising 
the wavelength of the Lyman-a break in the spectrum. 

Where redshifts have been determined by multiple 
groups, our results are generally in reasonable agreement 
with previous observations, given the difficulty in obtain- 
ing precise redshifts from the self-absorbed and resonantly- 
broadened Lyman-a lines at moderate spectral resolution. In 
the case of GLARE 1042, our measured redshift of z = 5.83 
is in close agreement wit h that measured by FORS 2 (z = 
5.82. 1 Vanzella et al.|l2004l ). by the GRAPES team {z = 5.8, 
Malhotra et alll2005l) and, using Keck/DEIMOS spectra, by 
Stanway et afl (|2004al . z = 5.83). GLARE 3001 was spectro- 
scopically identified by FORS 2 as a line emitter at z = 5.78, 
agreeing with our redshift of z — 5.79. 

GLARE 1054 and 103 were both placed at z = 5.7 
by GRAPES spectroscopy (|Malhotra et alll2005l ). This con- 
trasts with our measured redshifts of z = 5.9 for GLARE 
1054 and (tentatively) z = 7.4 for GLARE 1030. The dis- 
crepancy for GLARE 1054 is within the expected level of 
agreement for GRAPES grism spectroscopy and so these 
results are consistent. Given the low significance of our line 
candidate in GLARE 1030, and the i'-drop redshift selection 
function, the GRAPES redshift remains the more likely. 

The remaining source, GLARE 3000 was identified by 
the VLT/FORS2 observations as a Galactic star on the 
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basis of weakly detected [OI] and [Nil] lines. This source, 
w hich corresponds to t he unresolved candidate SBM03#5 
m IStanwav et aD l|2003h . falls in our 'possible' category. We 
note that the FORS II spectrum is flagged with their quality 
class 'C and that there is no secure continuum detection. 
Despite this, a Galactic star is still the most likely identifi- 
cation of this target, illustrating the caution with which we 
present our fainter line candidates. 

No other target on our 2004 GLARE mask has been re- 
ported as a line emitter by other teams, or has an estimated 
redshift from GRAPES spectroscopy. 



3.5 Limits on NV and other Emission Lines 

Lyman break galaxies at z ~ 3 show few other emission 
lines in the rest frame ultraviolet. Th e composite sp e ctrum 
of 1 000 such galaxies produced by IShaplev et al.l (120031 ) 
shows weak emission features due to Sill* (A rC st =1265A, 
I309A, 1533A), OIII] (A = 1661 A, 1666A) and CIII] 
(A =1908A), and absorption features due to stellar winds, 
primarily Sill and CIV. 

The presence of an AGN in our target galaxies could 
also lead to the presence of emission lines due to high exci- 
tation states, primarily NV (A rcst =1238. 8A, 1242. 8 A) and 
SilV (A =1394A, 1403A), and rarely OV] (A =1218A). 

For a galaxy at z = 6, our spectra extend to rest frame 
wavelengths of «1428A. Given that we are not able to mea- 
sure a high signal to noise continuum on any one target, we 
are unable to measure absorption features in the spectra, 
and therefore confine ourselves to searching for evidence of 
other emission features in the spectra of our Lyman-a emit- 
ters. 

A careful inspection of our five good Lyman-a emis- 
sion line candidates does not provide evidence for any other 
emission lines at the redshift of Lyman-a. While this is not 
surprising given the weakness of secondary lines, the failure 
to detect NV in this spectrscopy suggests a large Ly-a / NV 
ratio. Using our 3a limit on undetected lines as an upper 
constraint on NV flux we determine that /(Lya)//(NV) > 
[10.5, 4.5, 2.8, 5.1, 7.5] (3c) respectively for GLARE targets 
[1042, 1054, 1008, 3001, 3011] Typical line ratios for AGN 
are /(Lya)//(NV) = 4.0 (Osterbrock 1989), while those for 
galaxies are greater than this. The limits we determine dis- 
favour an AGN origin to the Lyman-a flux in all but the 
faintest target (in which the constraint is too weak to make 
a firm statement). These limits constrain AGN activity, if 
present, to only a weak contribution to the rest frame ultra- 
violet flux. 



4 THE EQUIVALENT WIDTH DISTRIBUTION 
OF GLARE LINE EMITTERS 

4.1 The Observed Equivalent Width Distribution 

Using the mean variance in the background of the exposed 
slits, and the broadband magnitudes of the targeted galax- 
ies, we are able to calculate limits on the rest frame equiva- 
lent width Wo for those sources which satisfy our colour se- 
lection criteria and yet are undetected in our spectroscopy. 
These are presented in table [5] In each case we assume that 
the galaxy lies at the mean i'-drop redshift (z — 6.0) and 



that the emission line has not been lost behind a the sub- 
traction residual of a bright night sky line. 

Combining our line limits with those sources for which 
emission lines have been identified or tentatively proposed 
forms a sample of 24 sources uniformly selected from their i'- 
drop colours. The resultant distribution of Lyman-a equiv- 
alent widths is plotted in figure [5] 

If all the sources for which candidate emission lines 
are identified in this paper prove to be z w 6 Lyman-a 
emitters, then the escape fraction of Lyman-a photons at 
z ~ 6 appears qualitatively similar to that at z ~ 3. From 
our z ~ 6 sample, 66% of the sources (16 out of 24) have 
Lyman- a equivalent w i dths < 25A, compared with 75% at 
Z Ri3 jShaplev et all l2003h . although the lower redshift 
sample loses a smaller fraction due to skyline contamina- 
tion. The samples probe to similar points on the luminosity 
function in both cases (approximately 0.1 L*). These frac- 
tions are consistent within the errors on our small number 
statistics. Harder to explain in comparison with lower red- 
shift galaxies is the tail stretching to very high equivalent 
widths (> 200A) observed in this survey, a trait also ob- 
serv ed in some narrowband sel ected sources at this redshift 
(e.g. IMalhotra fc Rhoaddliool) and in other ^-dropo ut Ly- 



man break galaxies (e.g. Dow-Hygelund et alj [2006, who 
find one source with Wo =150A). Although the number 
statistics are small, we observe four line emitters (17 ± 8%) 
with equivalent widths in the range 50 — 100A, and a further 
four with Wo >100A (th ree of which come from our 'robust' 
list of line emitters) . This contrasts with the Lyman Break 
Galaxy population at z ~ 3 in which < 5% of galax ies have 
line emission with Wo > 100A (|Shaplev et al.ll2003h . 

While high redshift galaxies at both redshifts are se- 
lected on their rest-frame ultraviolet continuum and the 
spectral break caused by Lyman-a absorption, the two pop- 
ulations are not identical. 

The sample discussed here reaches some two magni- 
tudes fainter than the tentatively pro posed and still uncer- 
tain typical lumino sity L* at z ~ 6 l|Bouwens et al.ll2006l ; 
iBunker et al.1 [20041 ). This contrasts wi th a sample reachin 
just one magnitude below L* at z = 3. Sh aplev et al.l (|200; 
considered subsamples at z = 3, dividing their spectroscopic 
data into quartiles based on rest-frame equivalent width. 
They found a modest trend in the strength of line emission 
with magnitude. Galaxies in their highest equivalent width 
quartile are some 0.2 magnitudes fainter on average than 
their quartile of weak to moderate line emission. It is possi- 
ble that the strong line emission observed here is more typ- 
ical of sources with faint rest-frame ultraviolet continuum 
suggesting that faint sources differ physically from brighter 
members of the population. 

A second difference between the samples is intrin- 
sic rather th an arising from a s el ection effect. Severa l 
authors fe.g. [ Stanway et al.1 120051 ; IBouwens et al.1 120061 ; 
ILehnert fc Bremerl 120031 7 have now observed that Lyman 
break galaxies at z > 5 ha ve steeper rest-f r ame ultraviolet 
slopes than those at z ~ 3. IBouwens et al.l (|2006l ) interpret 
this as indicating that the dust properties of this population 
evolve over redshift. While a steep rest-ultraviolet slope can 
also arise due to low metallicity or a top-heavy initial mass 
function (as discussed below) dust evolution is a natural in- 
terpretation. At z w 6 the universe is less than 1 gigayear 
old and galaxies may not have had time to develop a high 
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dust content. IShaplev et all (|2003h found that z ~ 3 galaxies 
with high equivalent widths in Lyman-a also had lower mean 
dust extinction. Lyman-a photons are resonantly scattered 
by dust and hence the line is preferentially absorbed with 
respect to the rest-frame ultraviolet continuum. If the distri- 
bution of Lyman-a equivalent widths in the z tn 3 popula- 
tion is truncated by dust absorption, this could produce an 
apparent 'excess' of strong lines at high redshift. However, 
even zero dust absorption cannot explain equivalent widths 
exceeding 200A unless the sources are also very young and 
very low in metallicity. The full explanation for the equiv- 
alent width distribution observed in the GLARE data may 
well be a combination of these effects and those discussed 
below. 

Several other possible explanations exist for both the 
steepening of the rest frame UV slope and the difference in 
equivalent width distributions. An interpretation of contam- 
inant galaxies at lower redshifts seems unlikely due to our 
i'-drop selection criteria; low redshift sources with strong 
spectral breaks are likely to have more than one emission 
line in our observed redshift range. 

A high equivalent width line can arise if the observed 
Lyman-a photons excited by a population of AGN rather 
than by young, hot stars. The luminosity function of AGN 
is poorly constrained at these magnitudes and redshifts, but 
the space density of such sources is expected to be very low 
(e.g., the z > 6 SDSS QSOs, Fan et al. 2003). At z « 6, the 
deep 2 Ms X-ray exposure o f the UDF and surrounding re- 
gion l|Alexander et al]|2003T ) would detect only the brightest 
AGN (L > 10 42 ergs s" 1 cm -2 ). None of the GLARE targets 
are detected in this X-ray observation. AGN would also be 
expected to show emission lines that are not just strong but 
also broad, while none of the GLARE line candidates are 
broad. There is also an absence of high-ionization lines such 
as NV 1240 A which are common in AGN. 

The tail of line emitters extending to higher equivalent 
widths may also arise from a difference in the fundamental 
properties of the stellar population between z w 6 and that 
at lower redshifts. 

Modeling of emission from metal-free Population III 
galaxies predicts rest frame Lyman-a equivalent widths 

> 1000 A for young starburs ts (<2 Myr), de creasing to W(Ly- 
a)~500A for older bursts fSchacrcr 200^). These very high 
equivalent widths arise from the relatively hard spectrum 
of metal-free reactions in the most massive stars. How- 
ever it is unlikely that zero metallicity (population III) 
stars are still contributing significantly to the emission from 
massive stars almost a billion years after the Big Bang, 
particularly given the identification of stellar populations 

> 100 Myr old in some z > 5 sources fe.g. lEgami et al . 2005; 
lEvles et alj|2005h . Further evidence for moderate metallic- 
ity at z > 5 has been observed in the spectroscopy of bright 
AGN from the Sloan D igital Sky Sur vey. Metals including 
iron l|Barth et al.l 120031 ') and carbon (|Maiolino et al.l 120051 ) 
have been detected from even the highest redshift quasar 
(at z = 6.4). 

By contrast, even l/20th solar metallicity leads to a 
sharp reduction in the peak (zero age) equivalent width pre- 
dicted to ~300A, w ith a more typical W(Ly-a )~100A by an 
age of 10-100 Myr (|Malhotra fe Rhoad j|2003 ) . Most of the 
candidate emission lines presented in this study can hence 
be explained with normal, if low metallicity, populations. 
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Figure 5. The distribution of rest frame equivalent widths for the 
GLARE mask, divided into bins of 25A. Line emission candidates 
are shown with solid lines, sources with equivalent widths > 150A 
are shown at 150A, and those for which only upper and lower 
limits are available are added to the distribution to produce the 
dotted lines. 



However at least two of our good line emission candidates 
have rest frame equivalent widths exceeding 200 A. This is 
possible if the galaxy is in the first few Myrs of an ongoing 
starburst, but may also provide evidence for variation in the 
initial mass function of star formation. 

High equivalent widths of Lyman-a emission can arise 
from a "top-heavy" initial mass function (i.e. star forma- 
tion weighted towards a p opulation of high mass stars). 
iMalhotra fc Rhoadsl (|2004bh calculated the Lyman-a equiv- 
alent width expected for a metal-enriched population with 
a very extreme IMF, weighted towards massive stars (i.e. 
IMF slope a = 0.5). As is the case for low metallicity popu- 
lations, the flux is weighted towards a harder spectrum, and 
Lyman-a emission is strengthened. They found that such an 
IMF could explain line equivalent widths of up to 240A at 
ages of a few Myr, with higher IMFs possible for very young 
bursts. 

The hard rest-frame ultraviolet spectrum associated 
with such an IMF may also be consistent with the steep rest 
frame u ltraviolet spectral sl opes observed in z ~ 6 galax- 
ies (e.g. IStanwav et ai1l2005h . While the evidence from the 
GLARE study is limited, with the number of high equivalent 
width sources small, the existence of such sources suggests 
that the environment of star formation at z ~ 6 is less metal 
enhanced, or weighted towards more massive stars than that 
at z ~ 3. 



4.2 Implications of High Equivalent Widths 

The presence of a tail of large equivalent width emission 
line sources has implications for the redshift distribution ex- 
pected from an i'-drop selection. 

The redshift distribution and number densities of i'- 
drop galaxies have usually been calculated assuming a neg- 
ligible contribution from line emission to the measured mag- 
nitudes and colours. In the case of line emitters with rest 
frame equivalent widths Wo < 30A this is a reasonable as- 
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ID 


Alt 


RA & Declination 


AB 


i'-z' 


W£f Q /A 


1001 


48989 


03 32 41.43 


-27 46 01.2 


28.26±0.12 


>2.1 (2a) 


<6.5 


1004 


46223 


03 32 39.86 


-27 46 19.1 


28.03±0.10 


>2.37 (2 a) 


<5.2 


1009 


12988 


03 32 38.50 


-27 48 57.8 


28.11±0.11 


>2.29 (2 a) 


<5.6 


1045 


21111 


03 32 42.60 


-27 48 08.8 


28.02±0.11 


1.67±0.26 


<5.2 


1047 


35271 


03 32 38.79 


-27 47 10.8 


28.44±0.14 


1.33±0.30 


<7.6 


1060 


11370 


03 32 40.06 


-27 49 07.5 


28.13±0.11 


>2.28 (2 a) 


<5.7 


1077 


16258 


03 32 36.44 


-27 48 34.2 


27.64±0.07 


1.42±0.16 


<3.7 


3002 




03 32 43.35 


-27 49 20.4 


26.89±0.07 


1.42±0.20 


<1.8 


3030 




03 32 48.94 


-27 46 51.4 


27.04±0.08 


1.41±0.23 


<2.1 


3033 




03 32 49.08 


-27 46 27.7 


27.18±0.09 


1.40±0.27 


<2.4 



Table 5. The 2004 GLARE targets for which no line emission was observed. Limits on the equivalent width are calculated using the z'- 
band magnitude to determine the continuum level, the 3 a standard deviation in the bac kground as the minim um line flux and accounting 
for IGM absorption, assuming the galaxy lies at z=6. 'Alt' indicates an alternate ID in lBunker et al. (2001) 



sumption, with the contribution to z flux amounting to a 
few percent. 

For line emitters with larger equivalent widths there can 
be a significant effect on the selection function of i'-drop 
galaxies. This effect falls into three main regimes as illus- 
trated by figure [6] If the emission line falls in the i'-band 
filter (i.e. z < 6.0), then the i — z' colour of the galaxy is re- 
duced by the line emission, and sources with high equivalent 
widths fall out of the colour-magnitude selection window. A 
source at our z detection limit with a rest frame Lyman-a 
equivalent width of lOOAcan be as blue as %' — z — 0.67 (AB) 
and an intrinsic line width of 150A would lead to a colour of 
just i' — z — 0.46 (AB). These colours are similar to those of 
much lower redshift galaxies. Ensuring a complete selection 
of line emitting galaxies at 5.6 < z < 6.0 is therefore im- 
possible using a simple colour-selected sample without also 
including a great many lower redshift contaminant sources. 

If the line emission falls in z' filter, the i'—z' colour is en- 
hanced and galaxies with continuum flux below the z' limit 
are promoted into the selection window. Hence at z > 6.5 
the i'-drop criterion can select a population of low contin- 
uum, strong line emission sources rather than the contin- 
uum break sources it targets. There is also a redshift range 
in which the line emission would fall in the overlap region 
between filters and both effects compete. 

This effect may lead to a bimodal redshift distribution 
for i'-drop galaxies, with weak emission line sources prefer- 
entially selected towards lower redshift, and high equivalent 
width sources selected at higher redshifts. It is necessary to 
quantify the equivalent width distribution, combining galax- 
ies at the faint limits explored here with brighter galaxies, 
in order to tightly constrain the galaxy luminosity function 
at any given redshift. 

As figure[6]illustrates, our good line emission candidates 
all lie in regions of parameter space affected by equivalent 
width selection biases. Our candidates lie in parameter space 
theoretically allowed for their equivalent width with two ex- 
ceptions - GLARE 1054 and 3011 both have strong emission 
lines in the tail of the i' filter that should have forced their 
i'—z' colours out of our selection function. We note, however, 
that that both of these sources have i' — z' colours within 
1 a of our minimum selection cutoff, and may have scattered 
up into out selection. Only one source is boosted into the 
selection function by virtue of its line emission rather than 
continuum flux, suggesting that our sample is, as expected, 
dominated by continuum-selected sources. If we are indeed 
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Figure 6. The effects of Lyman-a equivalent width on the i- 
drop selection function. Contours show the equivalent width line 
required to place a galaxy with the given redshift and continuum 
magnitude into the selection window. Solid lines indicate mini- 
mum equivalent width required, dotted lines indicate maximum 
equivalent width permitted. In the absence of line emission, all 
galaxies above the limiting observed magnitude and at z > 5.6 
are detected. In the presence of strong line emission, galaxies at 
z < 6.0 fall out of the selection as the emission line lies in the i'- 
band, while fainter galaxies at high redshift can be promoted into 
the selection due to brightening of the measured z' magnitude. 
The properties of the seven i'-drop sources with redshifts in this 
range are shown with diamonds. 



missing sources with high emission line equivalent widths 
at the low redshift end of our survey (where our redshift 
selection function is at its peak), then our conclusion that 
the distribution of line emission at high redshift extends to 
larger equivalent widths is strengthened. 

This effect has interesting implications for the luminos- 
ity functions presented in the literature for z ~ 6 galaxies. 
If the i'-drop criteria are overestimating the number density 
of faint continuum sources (due to contamination by a tail 
in the distribution of strong line emitters) then it is pos- 
sible that the faint end slope of the luminosity function is 
shallower than hitherto reported. If there is a large number 
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of strong emitters at 5.6 < z < 6.0 that have fallen out of 
the colour selection, then the number density of sources at 
these redshifts could be underreported. This may contribute 
towards the discrepancy between the observed evolution in 
the LBG luminosity function between z = 4 and z ~ 6 (e.g. 
iBunker et ail 12004 iBouwens et al.l 120041 ) and the reported 
lack of evolution in the Lyman-a emi tter luminosity func- 
tion o ver the same redshift range (e.g. Ma lhotra fc Rhoadsl 
l2004all . 



251) 



200 



Ando et al l (|2006l ) studied the equivalent width distr 



bution of a heterogeneous sample of i'-drop and narrowband 
selected galaxies with spectroscopically confirmed Lyman- 
a emission at z ~ 6. They commented on an apparent 
dearth of UVTuminous galaxies with strong Lyman-a emis- 
sion lines, noting an increase in the fraction of strong line 
emitters at Mab ~ —21. In figure[7]we show the equivalent 
distribution of the GLARE line emitters as a function of 
absolute magnitude. 

We detect five bright sources with a line luminosity 
brighter than 2 x 10 43 ergss -1 , thus occupying a regime un- 
occupied by the Ando et al sample. Of these, two possible 
emission line sources have equivalent widths greater than 
any in the previous sample. The four sources with a calcu- 
lated line luminosity between 5 x 10 43 and 10 44 ergss -1 are 
uniformly distributed in absolute magnitude. 

Although the number statistics of our sample is small, 
our sample is uniformly selected by broadband magnitude, 
removing possible biases due to the combination of two 
methods. Narrow-band surveys are biased towards sources 
with faint continuum and bright Lyman-a emission lines, 
whereas the i'-drop selection method is more uniformly sen- 
sitive to Lyman-a emission except at the very faint and low 
redshift ends of the sample. Our results do not support those 
of Ando et al., although clearly further observations are re- 
quired to clarify the fraction of line emitters at bright mag- 
nitudes. It is possible that the fraction of sources with high 
equivalent width line emission is subject to environmental 
effects (e.g. the local number density of galaxies). The small 
HUD F field is, of course, subject to the effects of cosmic vari- 
ance l|Bunker et al.l[20 04) . For sources as clustered as Lyman 
break galaxies at moderate redshifts are known to be, a vari- 
ance in number density of 40% is expected in a field this size 
jSomerville et alj|2004 ). 



5 DISCUSSION 

5.1 The Redshift Distribution of i'-drop Galaxies 

The redshift distribution expected for line emitters detected 
by the GLARE survey is the convolution of the selection 
functions imposed by the i'-drop selection criteria and by the 
spectral response of Gemini/GMOS in this configuration, 
modified by the probability of detecting a line that falls on 
a night sky emission line. 

In figure[8]we show the expected distribution of emission 
line redshifts, neglecting the influence of strong line emission 
on source magnitude and colour (see above), and assum- 
ing a Schecter (1976) luminosity function with parameters 
L* — L*(z = 3), a = —1.9 and </>* — 4>*{ z = 3) (Stanway 
2004). Changing the parameters of the luminosity function 
has a slight effect on the shape of the distribution, with a 
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Figure 7. The equivalent width distribution of GLARE line emit- 
ters as a function of absolute magnitude. The region of parameter 
space surveyed by the compilation of Ando et al (2006) is bordered 
by a dotted line. Dashed lines indicate the loci z = 6 emission 
lines with luminosities of 10 44 , 5 X 10 43 , 2 X 10 43 , 10 43 , 5 X 10 42 
and 10 42 ergss -1 (from top left to bottom right). Our equivalent 
width limit for sources without line detection corresponds to a 
rest frame luminosity of 10 42 ergss -1 , assuming those sources lie 
at z sa 6. 



reduction in L* with increasing redshift slightly broadening 
the peak of the distribution in redshift space. 

We take the conservative assumption that a source has 
a detection probability of 0% if lying on a very strong line 
(> 5x the inter- line sky continuum) and 50% on top of a 
weaker sky line (> 1.5 x sky continuum level). The night sky 
spectrum is measured directly from our spectroscopy. 

Clearly, our survey is most sensitive to sources lying 
towards the lower-redshift end of our redshift range, and to 
sources lying in the low OH line regions commonly exploited 
by Lyman-a surveys at z ~ 5.7 and z ~ 6.6. Nonetheless, we 
have a reasonable probability of detecting sources in the sky- 
line complexes at z < 6.5, particularly if the source emission 
lies between skylines. 

Given this theoretical sensitivity distribution, it is inter- 
esting to note that the five robust line emitters presented in 
this paper (and one possible emission line also lying in this 
redshift range) do not fall at the peak of the detection sensi- 
tiv ity distribution , but ra ther within the skyline complexes. 
In IStanwav et alj l|2004bl ) we suggested weak evidence for 
large scale structure in the HUDF field at z = 5.8, based 
on the three redsh ifts known at the time. More recently 
iPirzkal et al.l (|2004l ) have found evidence for an over-density 
at z = 5.9 ± 0.2 in the same field based on GRAPES Grism 
spectroscopy. Given that it unlikely that our fainter targets 
will yield redshifts with the current generation of spectro- 
graph, confirmation of this may prove challenging. However, 
our results so far appear consistent with large scale structure 
at z ~ 5.8 in this field. 



5.2 The Space Density of Lyman-a Emitters 

We have identified a total of 12 (5 of which were strong) 
candidate emission lines with a luminosity in the range 
(1 — 50) x 10 42 ergss -1 in our sample of 24 i'-drop galaxy 
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Figure 8. The redshift distribution expected for sources with 
line detection in the GLARE spectroscopy. Assumed a distribu- 
tion of magnitudes given by a Schecter luminosity function, with 
constant number density over the redshift interval 5.5 < z < 7.0, 
modified by the sensitivity of Gemini/GMOS in our configura- 
tion, and assuming that a line is undetected if under a skyline, 
and has a 50% probability of being detected if under a weaker 
skyline. The redshifts of our candidate line emitters are shown by 
asterisks, see section |5T] for details. 



candidates. If this is indicative of the i'-drop selected galaxy 
population as a whole, and if the equivalent width distribu- 
tion is invariant with magnitude, we would expect to ob- 
serve line emission of this strength in approximately 25 (11 
strong) of the 51 i'- drop galaxies identified in the HUDF by 
iBunker et all (|2004r ) (which were used by GLARE for candi- 
date selection). Hence, this suggests a line emitter space den- 
sity of 6.7(2.9) x 1(T 4 (AlogL)- 1 Mpc" 3 ), using the effective 
comoving volume for the i'-drop selection in the HUDF sur- 
vey of 2.6 x 10 4 Mpc 3 . This space density comparable to the 
space density of 6 .3 x 1 0~ 4 (AlogL) -1 Mpc -3 estimated by 
IShimasaku et all (|2006l ) for Lyman-alpha emitters of com- 
parable luminosity ((3 — 25) x 10 42 ) at z ~ 5.7 in the Subaru 
Deep Field. 

Similarly, at z = 5.7 the Lyman-Alpha Large 
Area Survey (LALA) found a number density of 3.2 x 
10~ 4 Mpc~ 3 for Lyman-cv e mitters with L > 10 42 ergs _1 
( Mal hotra fc Rhoadsl 12004 bl ). This would predict the detec- 
tion of 7 Lyman-a emitters to our detection limit in our sur- 
vey volume. This is consistent with, although slightly lower, 
than the number of emitters we detect. 

Inferring the space density of emitters from the rela- 
tively small sample in this survey is, of course, subject to 
large errors due to Poisson noise on the small number statis- 
tics. However, at least two physical explanations are possible 
for a discrepency between the line emitter density found by 
GLARE and by LAE surveys. As mentioned in section [STTl 
the HUDF is in a region with a known overdensity of sources 
at 5.7 < z < 5.9. GLARE is also likely to detect a higher 
apparent density of line emitters if the Lyman-a photon es- 
cape fraction is higher in the relatively faint galaxies studied 
here than in more massive galaxies. If this is true then the 
number of sources observed at faint magnitudes exceeds the 
number predicted from brighter surveys. This phenomenon 
may well be expected if the deeper potential well of more 
massive galaxies traps more dust and a denser ISM than 



that existing in the shallow potential wells of small, faint 
galaxies. 



5.3 The Emission line properties of sub-L* 
Galaxies 

This spectroscopy is the longest integrated exposure on a 
single mask using an 8m telescope. It reaches the faintest 
line limits available for a large sample with a homogeneous 
i'-drop colour cut. The 3 sigma limit in most slits is approx 
4.1 x 10" 19 ergcm- 2 s~ 1 A -1 . 

Hence we are able to observe line emissions from can- 
didates significantly fainter than other surveys targeting i'- 
drop galaxies. 

Given the expenditure in telescope resources required 
to reach these faint limits, it is likely that many of the 
HUDF candidates may well remain inaccessible to ground 
based spectroscopy until the advent of the Extremely Large 
Telescopes, and from space by JWST with NIRSpec. HUDF 
i'-drops can be identified with reasonable confidence to 
z' AB < 28.5 (Bunker et al. 2004). A source with z' AB = 28.5 
and an emission line of W=20A (rest) at z = 6, with com- 
plete Lyman-a absorption shortward of the line, would have 
a contamination fraction of 13.5% to the z'-band from this 
line, and a continuum flux of 6.4 x 10 -21 erg cm -2 s _1 A -1 - 
i.e. a continuum flux a factor of ~ 10 less than we are able 
to detect at the 1 a level per spectral resolution element 
(6.5 A) . However, such a source would possess a line flux of 
6 x 10 -18 erg cm -2 s _1 , well within our detection limits for 
emission lines. 

Hence, although emission lines are detectable in sur- 
veys such as this, even for the faintest of our target sources, 
continuum spectra will remain inaccessible for some time to 
come. Even GLARE 1042, the brightest z ft! 6 in the HUDF 
does not have sufficient signal to noise to in this spectroscopy 
to detect absorption features, even though the continuum is 
now detected at high significance. The detection of stellar 
absorption lines would provide valuable kinematic informa- 
tion on the galaxies and their outflows, however while such 
investigations may be possible on a stacked image of our 
faint spectra, it is impossible on individual slits. A stack 
of the faint GLARE spectra will be considered in a future 
paper. 



6 CONCLUSIONS 

In this paper we have presented spectroscopy of the faintest 
known sample of i'-drop galaxies, derived from the Hubble 
Ultra Deep Field. An exposure of 36 hours per target was 
obtained with the Gemini-South telescope. 

Our main conclusions can be summarised as follows: 

i) We have obtained extremely deep spectroscopy for 
twenty-nine science targets, reaching 3cr line flux limits of 
2.5 x 10 -18 erg cm -2 s _1 - corresponding to equivalent width 
limits of < 6 A for the majority of our targets. 

ii) We identify five i'-drop galaxies as good Lyman-a 
emitter candidates, four as possible candidates, and a fur- 
ther five tentative emission lines which we consider unlikely 
to be real. 

iii) We have considered the observed equivalent width 
distribution of i'-drop galaxies in the HUDF, and identify a 
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tail of line emitters with very high equivalent widths which 
is not seen in the lower-redshift Lyman-break galaxies at 
z « 3 - 4. 

iv) Several possible explanations for this effect exist. 
These include a tendency towards stronger line emission 
in faint sources, extreme youth or low metallicity in the 
Lyman-break population at high redshift, or possibly a top- 
heavy initial mass function. 

v) At the low redshift end of our selection function 
(5.6 < z < 6.0), the i'-drop selection method will fail to 
select line emitters with high equivalent width due to line 
contamination producing blue colours. In contrast, at the 
high redshifts end (z > 6), continuum- faint line emitters 
may enter the selection function if the Lyman-a line is suffi- 
ciently luminous. This has implications for the redshift and 
continuum magnitude selection function of i'-drop galaxies. 

vi) This sample significantly increases the number of 
faint i'-drop galaxies with known redshifts, and may begin 
to bridge the divide between continuum and line selected 
galaxies at z ~ 6. 
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